Spinons in the strongly correlated copper oxide chains in SrCu02- 
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We have investigated the spin dynamics in the strongly correlated chain copper oxide SrCu02 for 
energies up to > 0.6 eV using inelastic neutron scattering. We observe an acoustic band of magnetic 
excitations which is well described by the " Muller-ansatz" for the two-spinon continuum in the 
S=l/2 antiferromagnetic Heisenberg spin chain. The lower boundary of the continuum extends up 
to ~ 360 meV, which corresponds to an exchange constant J = 226(12) meV. Our finding that 
an effective Heisenberg spin Hamlitonian adequately describes the spin sector of this ID electron 
system, even though its energy scale is comparable to that of charge excitations, provides compelling 
experimental evidence for spin-charge separation. 
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The unique properties of the one-dimensional (ID) 
electronic systems in copper based chain materials con- 
tinue to attract theoretical and experimental attention. 
Not only are they a test bed for understanding the 
unusual electronic properties of high-T c superconduc- 
tors, they also allow experimental access to fundamen- 
tal physical phenomena in one dimension, such as the 
ID Mott insulator, spin-charge separation, non-Fermi- 
liquid (Luttinger liquid) behavior, and Peierls instabili- 

ties HEiH 

The corner-sharing chain cuprate SrCu02 and its sis- 
ter material Sr2Cu03 are of particular interest. They 
feature co-planar Cu04 square plaquettes arranged so 
Cu 2+ chains extend diagonally through rs 180° Cu-O-Cu 
bonds. Although the hopping integral t for the Cu 2+ 3d 
electron in this geometry is large, on-cite Coulomb re- 
pulsion U, stabilizes a Mott insulating state (MI) A 
similar bond arrangement occurs in the high-T c cuprates, 
where the corner-sharing chains form a 2D square lattice 
of corner-sharing plaquettes. An intermediate situation 
that may approximate the magnetic effects of stripes in 
high-T c materials _|6j , is found in the homologous series 
Sr„Cu„+i02ri+i .7|, where spin-chains form (n+l)-leg 
ladders. A "one-leg ladder", i.e. a single chain, is the 
simplest representation of the whole universality class of 
odd-rung spin ladders @- 

Super-sonic spin excitations are central to under- 
standin g th e low-energy electronic properties of cuprates 
0- EJ El- Spin waves in the high-T c parent mate- 
rial La 2 Cu04 have recently been characterized compre- 
hensively by magnetic inelastic neutron scattering (INS) 
[l2|. A two-dimensional (2D) dispersion relation corre- 
sponding to a nearest-neighbor Heisenberg spin coupling 
J Fts 140 meV and a four-spin cyclic exchange J c w 38 
meV was observed. The substantial cyclic exchange is a 
spectacular consequence of the electron itineracy in the 
underlying 2D Hubbard Hamiltonian. The significantly 



larger exchange interactions in the chain copper oxides 
as compared to their planar high-T c relatives is one of 
the mysteries in the electronic structure of cuprates El • 
An accurate, direct measurement of J is vital for solving 
this problem, as values determined by different indirect 
techniques differ appreciably. The temperature depen- 
dence of the magnetic susceptibility up to 800 K suggests 
J = 181(17) meV and J = 190(17) meV in SrCu0 2 and 
Sr2CuC>3, respectively El> while mid- infrared absorp- 
tion data yield an rs 37% larger value J « 260 meV [l4| . 
Such a large exchange constant would imply that the 
spinon- velocity matches the Fermi- velocity. The question 
therefore naturally arises whether or not spin excitations 
couple to charge excitations in these materials. Is it at 
all appropriate to describe a spin excitation in terms of 
a simple Heisenberg spin Hamiltonian when there is no 
energetic separation between the spin and charge sectors 
0? Here we report a neutron scattering study of spin 
excitations in the chain copper oxide SrCu02 which ad- 
dresses these issues and provides an accurate value for J, 
as sought by theorists for some time 0, El . 

SrCu0 2 has a centered orthorhombic crystal struc- 
ture (space group Cmcm) with lattice parameters a = 
3.556(2)A, b = 16.27(4)A, c = 3.904(2)A. The corner- 
sharing Cu-0 chains run along the c axis and come in 
pairs, stacked along the b direction El- The coupling be- 
tween chains in a pair proceeds through ss 90° Cu-O-Cu 
bonds and is extremely weak and frustrated 0,11]. There- 
fore, SrCu02 is essentially a single-chain compound, sim- 
ilar to Sr2Cu03. In both materials the inter-chain cou- 
plings are very small. In SrCu02 short-range static spin 
correlations develop below Tn ~ 5 K, and weak modu- 
lation of spin fluctuations in the direction perpendicular 
to the chains only occurs for energies E < 2.5 meV |l5j. 

Measuring spin excitations at energies as high as ~ 0.5 
eV by INS only recently became a realistic possibility 
with the development of the MAPS spectrometer at ISIS. 
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A SrCu02 sample with mosaic r\ w 0.5° and m r* 3.9 g, 
previously studied in Ref. 0, was mounted in a He-filled 
Al can, and cooled to T = 12(2) K in a closed cycle refrig- 
erator. The sample was aligned with its (hOl) reciprocal 
lattice plane horizontal, and with the c-axis (chain di- 
rection) perpendicular to the incident neutron wavevec- 
tor ki. It was fully illuminated by the « 45 x 45 mm 2 
incident neutron beam. We index wave vector trans- 
fer Q in the Cmcm orthorhombic reciprocal lattice as 
{Qin Qk, Qi) and define the equivalent in-chain wavevec- 
tor for the unit lattice spacing through q = 2ttQi. The 
magnetic cross-section was normalized using incoherent 
nuclear scattering from a vanadium standard. 

Data were collected for incident neutron energies Ei = 
98, 240, 517 and 1003 meV and with the frequency of 
the Fermi chopper chosen to have fairly coarse energy 
resolution (FWHM ~ 5 — 10% of Ej over the energy win- 
dows in Fig.l) in order to maximize neutron flux. This 
is important because the magnetic intensity is very weak 
for high energy transfers, where a large Q is required to 
satisfy the energy and the momentum conservation, and 
an exponentially small Cu 2+ form factor suppresses the 
magnetic scattering cross-section. The data for different 
Ei shown in Fig. 1, (a)-(d) focus on different ranges of 
energy transfer. The onset of the magnetic scattering at 
a highly dispersive lower boundary, ~ Jsin(g), is clearly 
observed in Fig. 1, (b)-(d). For E t = 98 meV, Fig. 1, 
(d), the steep dispersion is completely unresolved due to 
the coarse Q-resolution imposed by the sample mosaic, 
AQ m 7] ■ fcj. As a result, the scattering looks like a rod, 
centered for all energies at Qi — 2n + 1, where n is an 
integer (i.e. q = n(2n + 1)). The splitting of the rod into 
two branches becomes evident at E ~ 200 meV in Fig. 1, 
(c), and is apparent in Fig. 1, (b), which shows the top of 
the dispersion at E ~ 360 meV. On the other hand, the 
data for Ei — 1003 meV in Fig. 1, (a) clearly show that 
the scattering continuum persists up to an upper bound- 
ary, which is consistent with the dispersion ~ Jsin(q/2). 

The simplest framework for understanding the es- 
sential electronic properties of the cuprate materials is 
provided by an effective single-band Hubbard model 
E El El El El E3 In the ID case, relevant for the 
chain cuprates, the Hamiltonian reads, 

n e = -tJ2 (4,* c i+i,° + H - c -) + U Yl n ^ n 3'i • M 

},<* 3 

Here er = f, j, ,Cj i<7 are the electron cre- 
ation/annihilation operators at site j along the chain, 
and t and U are the nearest-neighbor hopping integral 
and the on-site Coulomb interaction, respectively. In the 
undoped stoichiometric materials there is one electron 
per site, ie the band is half-filled. 

For very large U the electron hopping is suppressed 
and the lattice sites are effectively uncoupled. The sys- 




Momentum Q| (units of 2tt/c) 

FIG. 1: Color contour maps of the normalized scattering in- 
tensity projected on the (Qi,E) plane, measured in SrCu02 
for (a) Ei = 1003 meV, \Q k \ < 7; (b) E, = 517 meV, \Q k \ < 5; 
(c) E, = 240 meV, \Q k \ < 5; (d) E r = 98 meV, \Q k \ < 4. 
An energy- dependent, but Q-independent non-magnetic back- 
ground scattering, measured at Qi w 0, was subtracted. The 
corresponding resolution-corrected intensity calculated from 
the Miiller ansatz Q is shown on the right, (e)-(h). 

tern for any spacial dimension is a narrow-band, wide- 
gap Mott-Hubbard insulator. In this limit, 2t/U -C 1, an 
energy-scale separation occurs, where a charge gap ~ U 
separates a ~ 2t wide conduction band from the low- 
energy, spin part of the electronic spectrum with much 
smaller bandwidth ~ 2t ■ (2t/U). The latter is described 
by an effective Heisenberg spin- 1/2 Hamiltonian, 

W = J^ Si s i+1 , (2) 

j 

with exchange coupling J = At 2 /U{l + 0(2t/U)) Ell- 
Such a situation occurs in many strongly correlated 
charge-transfer insulators, eg in KCUF3 |2lJ. In spite of 
the fundamentally different physical origin, it essentially 
resembles a band insulator: the electrons are strongly 
confined and the low-energy spin- waves are effectively de- 
coupled from the high-energy charge excitations. 

In the corner-sharing copper oxides U/t is not as large 
as for example in KCUF3, and the correlation energy is 
comparable to the hopping kinetic energy. Such proxim- 
ity to the Mott metal-insulator transition is believed to 
be crucial for the appearance of high temperature super- 
conductivity upon doping 0, llOj . A significant electron 
itineracy in 2D has immediate consequences for the spin- 
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FIG. 2: Constant-i5 cuts of the measured scattering intensity 
of Fig. 1 which show the lower boundary of the spinon con- 
tinuum. The intensity in each panel is combined within the 
energy range shown in the upper right corner. Solid lines are 
the fits to the resolution-corrected MA cross-section The 
broken lines are the calculated dCP boundaries Dashed- 
dotted lines show the Q-independent non-magnetic scattering. 

wave spectrum. This was recently observed in La2CuC>4 
P3, where U « 2.2 eV, t w 0.30 eV, and the band of 
spin excitations extends to w 0.30 eV (ie the nearest 
neighbor exchange coupling is J « 140 meV), which is 
non-negligible compared to the charge gap. As a result, a 
simple bilinear Heisenberg spin Hamiltonian such as Eq. 
(0) is inadequate, requiring non-physical values of the 
superexchange constants. An excellent account of the 
data was however obtained, by including four-spin cyclic 
exchange, J c ~ 0.3 J, which appears in second-order per- 
turbation theory for the 2D Hubbard Hamiltonian. 

The energy separation between the spin and the charge 
sectors of the electronic excitation spectrum is essentially 
absent in chain cuprates with w 180° Cu-O-Cu bonds 
|{|. Already, it is clear from the data of Fig. 1 that 
the top of the lower bound of the spin excitation contin- 
uum is at ~ 360 meV, resulting in J of ~ 220 meV in 
the spin Hamiltonian J2J. This translates into the total 
bandwidth for triplet spin excitations of ttJ w 0.72 eV, 
rou ghly equal to the charge excitation gap A c « 0.75 eV 
[lil Il7j. If the spin sector of the ID Hubbard Hamil- 
tonian in such a situation can still be adequately 
described by the simple S=l/2 Heisenberg spin Hamilto- 
nian J2J) , the spin excitation spectrum must be almost en- 
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FIG. 3: Cuts of the measured scattering intensity of Fig. 1, 
(a), which probe the upper boundary of the spinon continuum. 
Curves and annotations are the same as in Figure 2. 

tirely composed of pairs of the free, spin-1/2, "fractional" 
elementary excitations called spinons. They form a con- 
tinuum bounded by the des Cloiseaux-Pearson (dCP) dis- 
persion relations 



-J|sin(g)| <e(q) <tt.J\ sin(g/2)| 



(3) 



Although the exact ground state of J5J was determined 
a long time ago |23j |. an exact expression for the two- 
spinon contribution to the dynamic spin structure factor 
S(q,e) was obtained only recently [24j. The expression 
differs only very slightly from the approximate, semi- 
empirical "Miiller-ansatz" (MA) expression [25j . 



SW(<?,e) = 



A e(e-e L (q))0(eu(q)-e) 



2tt 



(4) 



which is based on consistency with sum rules and numer- 
ical simulations. Here, 9(x) is a step function, £L,u(q) 
are the lower and the upper dCP boundaries of Eq. J3J , 
and A ~ 1 is a prefactor introduced in Ref. [25] which 
we refine in a fit. MA J2J is routinely used to describe 
the two-spinon scattering, e.g. in KCUF3 [2]]]. Although 
it disagrees with the exact result [24| by predicting a 
step- like singularity at the upper dCP boundary , this 
artefact often disappears upon convolution with the in- 
strument resolution function. We compare our data with 
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FIG. 4: Parameters (a) J and (b) ^4 refined by fitting cuts at 
different E, including those in Figs. 2 and 3. The symbols, 
D , ■ , o, •. correspond to cuts taken with neutron energies 98, 
240, 517 and 1003 meV, respectively. 

the S MA {q,£), weighted by ( 7 r ) 2 iv|-|§ ^(Q)| 2 as ap- 
propriate for the scattering cross-section, and corrected 
for MAPS resolution H3- Here ( 7 r ) 2 = 0.290 barn, N 
is the number of Cu 2+ ions fcj, kf are the incident and 
the scattered neutron wave vectors, g rs 2 is the Lande 
g-factor, and F(Q) is the magnetic form factor. 

Each dataset in Fig 1, (a)-(d) was divided into a series 
of constant-E 1 cuts, which were independently fitted to 
the MA cross-section described above to obtain a pair 
of values A, J for each cut. A representative set of cuts 
with the corresponding fits is shown in Figs. 2 and 3. 
We obtain reasonably good fits, with A and J scattered 
within < 20% of their average values J = 226(12) meV 
and A — 0.55(6). The values and the error bars obtained 
from the fits are plotted in Fig. 4. The average values 
and the errors for A and J were obtained by perform- 
ing the weighted average of the data shown in the figure 
(for J the Ei = 98 meV data were not included, as the 
dispersion of the spinon dCP boundary is unresolved) . 

We found it very important to use the anisotropic mag- 
netic form factor appropriate for the Cu 2+ spins in the 
3d x 2_ y 2 orbitals [28|. If the spherical form-factor corre- 
sponding to the jo term only 28] is used, then the values 
of A for constant-energy cuts centered on similar ener- 
gies, but taken from data sets with different Ei, differ by 
up to a factor of > 5. This is because the same triplet 
energy is measured at different Q, making the data sen- 
sitive to the real space distribution of the spin density. 

Our findings show that up to at least 0.6 eV spin dy- 
namics in SrCu02 is indeed well described by the MA for 
the two-spinon triplet continuum appropriate for a spin 
Hamiltonian (J2J. This is also apparent from the remark- 
able similarity of the contour maps of the measured inten- 
sity in Fig. 1, (a) - (d) to the corresponding resolution- 
corrected intensity calculated from the MA cross-section 
and shown in the right panel of Fig. 1, (e) - (h). The 



experimental value for A is noticeably smaller than the 
value A = 1.347 for which Eq. accounts for the total 
spectral weight of a spin- 1/2 degree of freedom per site 
[25j. At present it is unclear whether this discrepancy 
results from limitations in high-Q form factor extrapola- 
tions or is a subtle effect of charge fluctuations. 

In summary, we have provided a detailed map of spin 
excitations in the chain cuprate SrCu02 for energies up 
to > 0.6 cV. Apart from possibly an overall scale factor, 
the inelastic magnetic neutron scattering data are indis- 
tinguishable from those projected for a spin-1/2 Heisen- 
berg Hamiltonian (J2J) with J = 226(12) meV. Because 
spin and charge energy scales coincide in SrCu02, these 
data provide experimental support for spin-charge sepa- 
ration in one dimensional Mott Hubbard insulators. 
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